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Abstract 

•/^ ■ We depict and analyze a new effect for wavepackets falling freely under a barrier 

1^ , or well. The effect appears for wavepackets whose initial spread is smaller than the 

> ■ nf 1 

vo ' combination \ -r-^, between the gravitational length scale /„ = — —7:7 and the 

^ : _ V l^ol (2 m2 5)1/3 

f-*) . initial location of the packet zq. It consists of a diffractive structure that is generated 

! by the falling and spreading wavepacket and the waves reflected from the obstacle. 

O 
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Q ■ length is included. 



The effect is enhanced when the Gross-Pitaevskii interaction for positive scattering 



The theoretical analysis reproduces the essential features of the effect. Experiments 



pi 1' emanating from the findings are proposed. 
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1 Introduction 



In the past years we have described a new phenomenon called : Wavepacket diffraction 
in space and time. The phenomenon occurs in wavepacket potential scattering for 

the nonrelativistic Schrodinger equation and for the relativistic Dirac equation. The effect 
consists in the production of a multiple peak structure that travels in space and persists. 
This pattern was interpreted in terms of the interference between the incoming spreading 
wavepacket and the scattered wave. The patterns are produced by a time independent po- 
tential in the backward direction, in one dimension, and, at large angles, in three dimensions. 
The multiple-peak wave train exists for all packets, but, it does not decay only for packets 

[w 

that are initially thinner than . — , where w is a typical potential range or well width and q 

V Q 

is the incoming average packet momentum. For packets that do not obey this condition the 
peak structure eventually merges into a single peak. The effect appears also in forward and 
backward scattering of wave packets from slits, 

The experimental breakthrough of Bose-Einstein condensation in clusters of alkali atoms|^ 
0] Q lead to the reinvestigation of the influence of the earth's gravitational field on the devel- 
opment of a quantum system. 

Gravity is currently being advocated as a mean to allow the extraction of atoms from the 
condensate for the realization of an atom laser continuous output coupler . [|, |^, [1^ Despite 
the weakness of the gravitational force on earth, it has a major influence on atoms that are 
cooled to microKelvin temperatures. It is then necessary to include the effects of gravity in 
theoretical calculations with condensates. 

Many other gravitational effects with quantum systems are being considered nowadays. 



such as bound states of neutrons in a gravitational field above a mirror, [11], or the use 



of the coherence properties of condensates to serve as interferometers in the presence of 



gravity 12, 13 



The Bose-Einstein condensate in a magnetic trap is in reality a wave packet. To the 
extent that decoherence effects are not dominant, it is expected to evolve in a gravitational 
field in the same manner as a Schrodinger wave packet. In the present work we investigate 
the effects of gravity on falling packets. 

It will be shown numerically and analytically that packets falling under an obstacle but, 
free from below, display distinctive quantum features due to their wave nature. 

The Schrodinger wavepackets not only fall, in accordance with the equivalence principle, 
but also spread. The thinner the initial extent of the packet the broader the spectrum of 
momenta it carries. Consequently, it will generate many more components able to reflect 
from the obstacle, be it a well or a barrier. These reflected waves will interact with the 
spreading and falling packet. 

There is crossover length scale, at which the interfering pieces start to produce a diffrac- 



comprehensive bibliography on Bose-Einstein condensation may be found at the JILA site 
http://bcc01.phy.GaSoU.edu/bec.html/bibhography.htmj 
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tive coherent structure that travels in time, analogous to the effect of wave packet diffraction 
space and time previously mvestigated.i-i 

This length scale is the gravitational scale L = -, ^ — 777^. For Sodium atoms it is 

(2 fif) ' 

about 0.73 microns, while for Hydrogen it is 5.86 microns.0 For packets initially narrower 
than a proportion of this scale, the effect is extremely evident, and gets blurred the wider 
the initial packet. This effect may be observed with the same setup as the one used in 
Bose-Einstein condensation experiments, provided pencil-like, thin packets are produced 
and allowed to fall under a roof. [0 

We will provide analytical approximations to the exact solution of the problem that 
reproduce quite satisfactorily the numerical results. In section 2 we present numerical results. 
Section 3 will deal with theoretical aspects. Section 4 summarizes the paper and provides 
concluding remarks regarding possible experiments. 



2 Packets falling under a roof 



Matter wave diffraction phenomena in time induced by the sudden opening of a slit, 
or in space by fixed slits or gratings, are understood simply by resorting to plane wave 
monochromatic waves. 



Atomic wave diffraction experiments ||T6|, have confirmed the predictions of diffraction 
in time||14|| calculations. These patterns fade out as time progresses. 

The phenomenon of diffraction of wavepackets in space and time was presented in 
It consists of a multiple peak traveling structure generated by the scattering of initially thin 
packets from a time independent potential, a well, a barrier, or a grating. The condition for 
the pattern to persist was found to be 



a « - (1) 

V Qo 

where a is the initial spread of the packet, w is the width of the well or barrier and go is 
the impinging packet average momentum. For packets broader than this scale the diffraction 
pattern mingles into a single broad peak. 

The original motivation for the present work, was the addition of gravity to the poten- 
tial affecting the packet propagation. As described above, the effects of gravity become 
increasingly relevant to the dynamics of packets in traps and elsewhere. 

The educational literature abounds in works dealing with the dynamics of packets falling 
on a mirror. The so-called quantum bouncer is a clean example of the use of the Airy 



packet in the treatment of the problem. The use of the Airy packet is straightforward above 
the mirror with the boundary condition of a vanishing wave function at the location of 



^ We use fi. = 1, c = 1, and units of length in microns, of time in milliseconds, g = 9.8- — - 

(msec) 
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the mirror, and, becomes a nice laboratory for the investigation of the quantum classical 
correspondence, revivals, the Talbot effect, etc. Falling packets were not studied, perhaps in 
light of the preconception that nothing interesting will be found besides the expected spread 
and free fall of the packet. 

However, there is a surprise awaiting us here. This is not totally unexpected due to the 
wave nature of the packet that consists of modes propagating in both the downwards and 
the upwards direction. Analogous effects are apparent also when a packet propagates in 



parallel to a mirror without ever getting close to it.|]18| Again the spreading and interference 
between the incoming and reflected waves produce a wealth of phenomena. 

In this section, we present the numerical results for the falling of packets under a barrier 
or well using Gaussian wave packets. The use of Gaussian packets permits a straightforward 
connection to theoretical predictions. 

The scattering event starts at t =0 with a minimal uncertainty wavepacket 



^0 = ^ e'" 
a = qiz-zo) (2) 

centered at a location zq large enough for the packet to be almost entirely outside the 
range of the potential, a denotes the width parameter of the packet, q = mv is the average 
momentum of the packet. The potential affecting the packet is a square well, the gravitational 
interaction, and eventually the Gross-Pitaevskii (GP) interaction, that subsums the effects 



of forces between the atoms in the condensate in the mean field approximation [19 



V = mgz + U Q{w/2-\z-w/2\)+ ^ (3) 
where m is the mass of the atom taken here to be Sodium, w is the width of the well or 



barrier, of depth or height U, and g is the strength of the GP interaction |]19[, which using 
the scattering length of Sodium and a typical number of atoms in a trap of the order of 250 
atoms//i^ yields g ^ 25 for a wavefunction normalized to one. 

Figure 1 depits the setup of the problem. A thin wave packet is initially located around a 
point Zq = —7n under a thick barrier, the mirror. The mirror could be any fiat surface. We 
have represented the strength of the barrier by means of a square located above z = 0. The 
simplification of the present investigation assumes a the mirror to be a finite width plate of 
thickness t = 10/i in the figure extending to a much larger distance along the x,y plane. For 
the one dimensional calculation we took a thin packet along the z axis of width 4a ~ 0.4/i. 
It is supposed to be a pencil-like packet whose extension along the x,y plane is larger than 
the width , such that the quantum dispersion affects primordially the behavior of the packet 
along the vertical direction. In this direction, the facket falls freely. The interference between 
the packet components reflected from the barrrier, upwards moving waves, and those that 
fall freely without reflection, the downward moving ones, will induce a diffraction pattern. 
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Figure 1: Setup of the problem. The packet's initial location is = — 7/i. The mirror is 
depicted as a square of arbitrary height representing the strength of the repulsive potential. 
The mirror is located above z = and its width is 10 fi 

We will find that the contrast of the pattern depends crucially on the initial location of the 
packet and its width. 

An ideal implementation of this setup would be to use a Bose-Einstein condensate opti- 
cally trapped in a region below a plate. Suddenly stop the confining potential and let it fall 
freely under the plate. 

The condition for the observation of the effect would be to produce a thin enough con- 
densate. For repulsive Gross-Pitaevskii interactions as in Sodium, Rubidium, or Hydrogen, 
the width limitation is less stringent, as will be shown below. The extra repulsion of the 
interaction, on top of the quantum dispersion, effectively pushes the waves inside the con- 
densate against the wall. It acts as if the packet was initially thinner than it actually is. In 
particular, for Hydrogen the initial width of the packet may be as large a few microns. In 
section 4 we make some further comments concerning the relevance of the effect for atom 
lasers. 

The algorithm for the numerical integration of the Schrodinger equation of the present 
work is described in previous works. ^, Flux conservation for initially normalized 
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packets and energy conservation require that 



E 



-oo 
oo 



1 d'^{z,t)d'^*{z,t) 



2 m dz dz 2 



+ {mg z+ V{z) + ^ l^(^,i)r) l^(^,^)P 



dz (4) 



with E a constant independent of time. 

The numerical runs presented below achieved an accuracy in the flux conservation of 
around 0.2%, while the accuracy in the energy conservation was around 2%. The wave 
function was also found to obey the Schrodinger equation to an accuracy better than the 1% 
level even at points near the edges of the integration range. 

The only length scale appearing in the problem is easily derived from the Schrodinger 
equation to be Ig. For Sodium this is Ig = 0.73/i 

We consider the free fall of packets with initial widths a smaller and larger than a fraction 
of Ig with or without the GP interaction. We take the barrier to have the fixed strength 

4 TT (2 

oi U = lO^sec"^ obtained from U ^ = A^, with a, the scattering length of the 

m 

Sodium-solid scattering and N the density of a typical solid. For the well we use a value taken 
from a van der Waals type of strength at a distance of 1 nm, namely U ~ — lO^sec"^. 
We used a very high value for the attractive well strength, beyond the limit of applicability 
of the Lennard- Jones formula[^, to see whether even a large and unrealistic value for the 



attractive potentials influences the results as compared to the repulsive case. 

Figure 2 shows the time evolution of a thin packet falling under a mirror. This is essen- 
tially what is expected to be observed in an actual imaging of a Bose-Einstein condensate 
falling freely under a mirror as a function of time. The figure shows the results for Sodium, 
however, as mentioned above the same will occur with a Hydrogen condensate, with much 
wider packets initially. 

Figure 3 shows wide and thin packets profiles after 4 msec fall under a repulsive barrier. 

The behavior of a thin packet is qualitatively different. A wide packet falls undistorted 

except for the natural spreading. A thin packet whose width is smaller than J ^ (see below 

for the appearance of Zq in the expression), possesses a distinctive diffractive structure. The 
rightmost (upper) edge of the structure resembles the Airy packet absolute value, however, 
the packet drops exponentially at large \z\ values, whereas the Airy packet diminishes as 
In the next section we will address a theoretical approach to the problem. Approxi- 
mate analytical solutions will be provided that reproduce the basic features of both the thin 
and wide packets. 

In figure 4 we present an analogous picture for the case including the GP interaction. As 
expected ||19|| , the distinctive feature of the inclusion of this interaction is an extra repulsive 



force, due to the positive scattering length, that enhances the broadening of the packet, 
and consequently the diffraction effect. As seen from the figure, for the same values of 
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Figure 2: Time evolution of a packet profiles of initial location zq — — 7 falling under 
the repulsive barrier (not depicted). 
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Figure 3: Packets profiles of initial location zq — — 7 after t— 4 msec falling under a 
repulsive barrier. 
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Figure 4: Packets profiles of initial location zq — — 7 fj, after t= 4 msec falling under a 
repulsive barrier and subjected to the Gross-Pitaevskii interaction 



initial packet widths, the nonhnear repulsive interaction produces a much cleaner interference 
pattern than the case lacking it in figure 3. 

Figure 5 depicts the influence of the GP interaction on the packets profiles for initially 
thin packets. 

The GP force produces a diffractive structure that has a much starker contrast. Large 
momenta arc excited by the repulsive interaction, producing a more effective superposition 
between incoming and reflected waves. 

We argued that the appearance of a diffractive structure is determined by the ratio 
e = — Figure 6 shows the results when this ratio is around 1. While the pattern has 

almost disappeared from the faUing packet devoid of GP interaction, it has lost contrast, 
but not disappeared completely from the packet subjected to the GP force, e is then the 
relevant ratio for the appearance of the diffractive structure. 

Replacing the barrier by a well has no effect whatsoever for packets without initial mo- 
mentum. The higher the initial momentum of the packet (if positive), the easyer the trans- 
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Figure 5: Thin packets profiles after t= 4 msec of fall under a repulsive barrier with and 
without the Gross-Pitaevskii interaction. 
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Figure 6: Borderline packets profile after t — 4 msec of fall under a repulsive barrier without 
and including the Gross-Pitaevskii interaction. 
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mission through the welL The differences arise then for packets thrown against the obstacle 
only at large momenta as compared to This point will be dealt with in the future. 

3 Theoretical approach to faUing packets 



The Airy function |21] is the solution of the Schrodinger equation for the propagation in a 
uniform field. Despite being named Airy packet, it is not normalizable and belongs to the 
set of wave functions in the continuum. It is the analog of a plane wave in free space. 

The educational literature abounds in references to uses of this packet, both in the context 



of the quantum bouncer^T^ and in the treatment of generalized Galilean transformations. p2 
For the solution of a packet falling on a mirror, the quantum bouncer, the Airy functions 
shifted to positions that have a zero at z = 0, serve as a basis to find the time development 
of an arbitrary initial packet. The non-normalizability of the Airy packet is of no hindrance 
here, because only the upper decaying part of the packet is used. Note however that an 
aspect that is ignored in the literature is the absence of orthogonality between the different 
shifted packets when integrated only over the positive z axis. This is perhaps not a severe 
problem, but was not taken into account in the works using the Airy packet for the quantum 



bouncer problem. ||T7| 



The lower piece of the Airy packet is oscillatory and does not decay fast enough to serve 
as a basis for packets initially located under the mirror. Only when a continuum of energies 
(both positive and negative) is used, it is possible to expand an initial wavepacket located 
at z < in terms of Airy functions. 

Before analyzing the problem by means of Airy functions, we will develop a very simple 
approximation that allows the identification of the relevant length scale for the effect to 
occur. 

We will resort to a set of solutions that connects directly to plane waves. This solutions 
are not stationary in the rigorous sense of the word, because they have a time dependent 
phase that is not linear in time. The solutions are derived by transforming a plane wave to 
an accelerating frame p2|, namely 



Xk{z,t) = e'^ 



(f) = g z t+ k { z+ — ) - ^ ^ q — (5) 

Direct substitution in the Schrodinger equation proves that this family of solutions solves 
the equation for a potential V = m g z. 

Given that the initial wavepacket of eq.(|^) may be expanded readily in plane waves 
that coincide with eq.(^ at t=0, the Schrodinger equation then insures that the subsequent 
propagation of the packet will be obtained by replacing the plane waves by the solutions in 
the gravitational potential of eq.(|^). 

We find the expression at all times for a freely falling normalized gaussian packet 
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Figure 7: Profile of a packet of initial width a = 0.35 fi and initial location zq = — 7 
in free fall after t = 4 msec. Numerical calculation, dotted line, and theoretical formula of 
eq.(|^ , full line 
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Figure 7 depicts the comparison between the expression of eq.(|]) and a numerical calcu- 
lation. The agreement is satisfactory without any rescaling. It lends confidence in both the 
numerical scheme and the theoretical formulae. 

At t=0 we can readily build a packet that solves the problem with the boundary condition 
of $(0) = 0, corresponding to an impenetrable mirror at z = 0. Just an image packet located 
at —Zo in the inaccessible region above the mirror will serve. The solution at t=0 then 
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becomes ^{z,t = 0) = "^{z — Zo,0) — "^{z + Zo,0), with \I/ in eq.(||). $ obeys the equation 
of motion and the boundary condition. It also coincides with the initial packet of eq.(0) in 
the allowed region of z < 0. Propagating $ forward in time we obtain 



^z,t)^^iz-zo,t)-^iz + zo,t) (7) 

We write the approximate sign, because the cancellation at z = is only effective at 
short times. As soon as t increases, the wave function of eq.(|^ does not vanish any more. 
Presumably, an infinite set of image packets is needed. 

We could not find a closed analytical solution at all times. The solution of eq.(|^) repro- 
duces reasonably well the falling packets and distinguishes clearly between a packet narrower 
than Ig and one wider than Ig for zq of the order of a few packet widths. In order to see 
this we write the absolute value of eq.(0) for a packet with initial momentum g = for 
t » 2 m a"^ 



\^{z,t)\ = A e'^' ^ sin^{e2) + sinh^{e^) 

((. + g t^/2f + zi) 



02 



4t2 

m Zq {z + g 



III u yx. ^ y L I z.) 

^3 = 



t 

zo {z + g t^/2) 
YT^ 

63 is responsible for the blurring and loss of contrast of the diffraction pattern determined 

by the sin function. The larger 6^ the less visible are the oscillations. The criterion for the 

visibility of the pattern may then be written as max{63) << 1. The maximum value of - z 

g ^2 

is given by the descent of the packet and its spreading. Both are of the order of We 

can then write the condition for the visibility of the interference fringes to be 

,^ , a"^ Zq q ^ 
max(^3) ~ "^"^ ^ 

However, typically zq amounts to a few times the width of the packet, otherwise the 
oscillations will have a very large wavenumber, and will blur the pattern anyway. 
Hence we can write eq.(P) as 

//, N a"" g 
max\P^) ~ ^ << 1 

<< 1 



4 P 



or 

a « Ig (10) 
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If we keep zq, eq. (|TO|) becomes 



a « 2 J-S- (11) 

V ^0 

Eq.(pTD demonstrates that the relevant borderhne between a visible and blurred packet 

[W 

is a = \ —. A packet initially narrower that Ig located at a distance of a few times its 
V 

width under a mirror will definitely display interference fringes. Eq.(^ also tells us that this 
pattern travels with the packet unscathed. 

For long enough times we can improve upon the solution of eq.(|^) in order to compensate 
for the inaccuracy in the cancellation at z = 0. The correction is achieved by multiplying 
the subtracted wave by a space independent, but time dependent admixture factor, 



<^{z,t) ^ "^iz- zo,t) - X'${z + Zo,t) 

A = (12) 

By adding this factor we spoil the solution. Eq. ( [T^ ) does not solve exactly the Schrodinger 
equation, whereas eq.(|^) does. However, this inaccuracy decreases as a function of time, 
because A — 1 as t — ^ oo. For times not so long it improves a little the agreement with the 
numerical results. It apparently compensates for the need of an infinite set of packets that 
insure the boundary condition at all times. 

Figures 8 and 9 show the comparison between the formula of eq.(|12D and the numerical 
results both for a thin packet and wider one. 

The agreement is reasonable quite good for the wide packet and qualitatively reasonable 
for the thin packet. The formula captures the gross features, such a as the absence of a 
diffractive structure for a wide packet and the wavenumber of the oscillations for a thin 
packet. It has one evident limitation: Lack of a sharp cutoff of the packet at small distances 
-z for a thin packet. Nevertheless, the simple picture of single image packet is essentially 
correct. 
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Figure 8: Thin packet profile of widtli a = 0.3 /i after t = 4 msec falling under a reflecting 
mirror. Numerical calculation, solid line, and theoretical formula of eq. ( 1121 ) 
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Figure 9: Wide packet profile of width a = 2 fi after t = 4 msec falling under a reflecting 
mirror. Numerical calculation, solid line, and theoretical formula of eq. ) 
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We now proceed to find a solution in terms of Airy functions.^ 

A solution of the mirror problem for 2; < can be determined by means of a set of 
functions that solve the freely falling packet, namely linear combinations of the independent 
Airy functions Ai and Bi^^. A set that implements the boundary condition of a vanishing 
wave function at the origin reads 



Xa{x) = {Bi{a) Ai{x + a) - Ai{a) Bi{x + a))/^A{^{a) + Bi^{a) (13) 
z —E 

where x = —, a = — , with E, the energy of a stationary solution —00 < a < 00. 

Ig mglg 

The functions obey 



Xa{x) Xb{x) dx = 5{a-b) (14) 

-00 

with S, the Dirac S function. The set is orthonormal and complete. 

The initial wave packet of eq.(@) is expanded in terms of the set above as 

/oo 
C{a)xa{x) da 
-00 

C{a) = f Xa{x)^{x,t = Q) dx (15) 
J —00 

The wave function at all times then becomes 

/oo 
C{a)xa{x) e'^' da 
-co 

(16) 



For thin enough wave packets C{a) may be obtained analytically by extending the inte- 
ition to +00 and neglecting corrections of order 

packet. Under these approximations C(a) becomes 



gration to +00 and neglecting corrections of order 7 = with a, the initial width of the 



2 7^ 



C{a) = [2-'nYj Ig exp{{a + xo) Y + ^) Xia + xo + r) (17) 

£0 
I, 

Inserting this formula into eq.(0) we find the wave function at all times. 



where Xq 

^9 



•^The author wishes to express his gratitude for the extremely helpful remarks of one of the anonymous 
referees concerning the use of Airy functions 
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Figure 10: Numerical proGle, dashed line, and, analytical solution of eq.dl 
the thin packet of figure 1 
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Figures 10 and 11 present the comparison of the numerical calculation for a thin packet 
and a wider packet to the analytical formula of eq.(|T6|) with the approximation of eg. (p!7D . 
The thin packet fit is quite good, while for the wider packet there is a discrepancy in the 
height of the peak. This discrepancy is due to the fact that 7 for a wider packet is not 



negligible and the formula of eq.(17) needs to be ammended. However in both cases the 
basic feature of the existence and absence of a diffractive train is evident. It is worth 
mentioning that there is no rescaling of neither analytical nor numerical data points. 

The long time behavior of the thin packet as compared to a wide packet remains un- 
changed for as far as we could integrate the time dependent Schrodinger equation numeri- 
cally and still agrees with the analytical solution of eg . ([TBI) . We reached a time of 30 msec 
and the profiles just spread, but do not change in shape. At that time the center of the 
packet is around 4.5 mm below the mirror. 

The integration of the equation becomes prohibitive beyond this time due to the strong 
oscillations in the wave function that require an increasingly smaller time step. The results 
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Figure 11: Numerical profile, dashed line, and, analytical solution of eq.(]7^j, solid line for 
the wider packet of figure 1 



20 



however look quite firm. The diffractive structure persist to infinite times, as it was found 
for the wavepacket diffraction in space and time effect. [|T|-|^ 

In the next section we will provide some closing statements on the relevance of the falling 
packet effect for Bose-Einstein condensates and atom lasers. 



4 Summary 

We have found that a falling packet that is blocked by an obstacle from above, behaves 
differently depending on its initial spread. A packet wider than the gravitational mass 
dependent length scale Ig falls almost as a free packet except for a lag due to an effective 
attraction to the barrier or well, while a thin packet has a distinctive diffraction pattern that 
propagates with it analogous to the 'Wavepacket diffraction in space and time effect '|[T|]-[^. 

We consider now a possible scenario to implement the findings of this work, and also 
take advantage of them for the purposes of creating an atom laser. Although there seems to 



be some controversy as to what an atom laser is[^, at a pedestrian level it would consist 
in a three stage machine: A feeding stage that pumps in atoms in an incoherent phase; a 
condensation cavity and a continuous output coupler. The investigation of these stages is as 
of today very advanced both theoretically and experimentally. 



We would like to point out that the results of this paper suggest an alternative avenue 
for the continuous output coupling of a Bose-Einstein condensate. 

Basically, it consists of an orifice, a physical one or one drilled in the mesh of laser 
radiation that confines the condensate. Through the orifice the condensate can exit in a 
pencil-like thin jet of atoms. 

Position now a mirror above the atoms and let them fall freely while the feeding continues. 
It appears then that the outcome will be a coherent train of atoms having the characteristic 
oscillations found here, provided the width of the pencil is smaller than ^I^/zq. 

In order to see if this design works and check whether it fits the criteria of an atom laser 



set in ref.|27], we need to perform at least a two-dimensional calculation with a source term. 



This endeavor is currently underway. 
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